Bipartite bound entanglement in continuous variables through deGaussification 
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We introduce a class of bipartite entangled continuous variable states that are positive under partial transpo- 
sition operation, i.e., PPT bound entangled. These states are based on realistic preparation procedures in optical 
systems, being thus a feasible option to generate and observe genuinely bipartite bound entanglement in high 
precision experiments. One fundamental step in our scheme is to perform a non-Gaussian operation over a 
single-mode Gaussian state; this deGaussification procedure is achieved through a modified single-photon ad- 
dition, which is a procedure that has currently being investigated in diverse optical setups. Although dependent 
on a single-photon detection in a idler channel, the preparation can be made unconditional after a calibration of 
the apparatus. The detection and proof of bound entanglement is made by means of the Range Criterion, theory 
of Hankel operators and Gerschgorin Disk's perturbation theorems. 

PACS numbers: 03.67.Mn, 03.65.Ud, 42.5().Dv 



The characterization of the phenomenon known as bound 
entanglement [l] constitutes one of the greatest challenges 
in quantum information theory |2|. Bound entangled states 
are those that, although being entangled, do not allow dis- 
tillation of any pure entangled state with Local Operations 
and Classical Communication (LOCC). This kind of entan- 
glement is associated with non-intuitive theoretical aspects of 
quantum information processing, such as communication us- 
ing zero capacity channels |,3J, or the irreversibility of entan- 
glement under LOCC 10]^. They also have practical impli- 
cation in quantum cryptography fSl, channel discrimination 
[7 1 and many quantum information protocols in general |9|. 
Experimental realization of such states has been only recently 
achieved ||T0HT3]| and is restricted until now to the multipartite 
scenario. Particularly, for continuous variable systems, local 
Gaussian operations, which are relatively simple to implement 
experimentally, are not able to distill entanglement if the sys- 
tem state is Gaussian 1,14., 15 1. In fact it is impossible to gen- 
erate bound entanglement with bimodal Gaussian states lfT6l . 
Quite recently an experimental investigation explored this fact 
for the generation of such a kind of states employing a bipar- 
tition with more than two modes at Alice and Bob sides ifTSl . 
To generate a two-mode continuous variable bound entangled 
state one necessarily has to move it out from the Gaussian 
class of states, by implementing some non-Gaussian opera- 
tion over a Gaussian state (deGaussification). 

In this article we propose a class of genuinely (two-mode) 
bipartite bound entangled states of arbitrary dimension that 
can be unconditionally prepared in optical systems with sim- 
ple extension on current experimental techniques. Our ap- 
proach is to deGaussify a thermal state by a photon-addition 
followed by an incoherent mixing with a squeezed vacuum 
in an orthogonal polarization. These states serve as the in- 
puts for generation of entanglement when mixed with vac- 
uum states. The detection of entanglement is simple and is 
made using the well-known Range Criterion ifTTl . Choosing 
properly the parameters involved, it is possible to design a 



state that is Positive under Partial Transposition (PPT), hence 
undistillable [18|, representing a practical improvement over 
the findings of Ref. 1 19 1. As a side result we derive some the- 
oretical insights on one-mode states through connections be- 
tween PPT property and Hankel operator theory (see Propo- 
sition 1), which are used together with Hadamard products 
and Gerschgorin Disk's perturbation theorems for the proof 
of bound entangled states. 

Preparation. It is impossible to generate a bound entan- 
gled state for a system composed of a single mode for Al- 
ice and an arbitrary number of modes for Bob, if their joint 
state pAB is Gaussian - Positivity of the partial transpose of 
PAB sufficiently implies separability |fT6l . The sufficiency 
no longer holds if Alice has more than one mode, such as in 
the experimental implementation on Ref. 1 13|, or if the joint 
state is not Gaussian. The procedure employed here for gen- 
erating bound entanglement between two modes employs the 
later approach, and requires two fundamental steps - a photo- 
addition over a thermal state, and an incoherent mixture to a 
two-mode squeezed vacuum. The experimental setup is pre- 
sented in Fig. [T] There a V^-polarized single-mode is prepared 
in a thermal state, pT ~ X]ra?'"-I'^)("'l' ^^'■^ thermal distribu- 
tion {pn}, through a phase-randomization process by a rotat- 
ing disk (RD)|29|. This light mode. A, is then photon-added 
(through a process inside the box PA to be later described), 
transforming {p„} into {p^}, and then mixed with a vacuum 
state mode, B, on the 50 : 50 beam-splitter BS-1. Since an 
arbitrary beam-splitter action on the creation operators of two 
input modes is given by the global unitary operation 
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the output state for the 50 : 50 (6* = 7r/4) beam-spUtter BS-1 
is given by 
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FIG. 1 . (Color Online) Two orthogonally polarized components of a 
laser beam are split in the polarizing beam splitter PBS- 1. The V- 
polarized beam is prepared in a thermal state pr through the rotat- 
ing disc RD and photon-added in the process described in the box 
PA (see text for details) to originate state pi. This mode is sent 
through the 50:50 beam splitter BS-1 whose output is described by 
Eq. pi. The iJ-polarized component of the laser beam coming out 
of PBS- 1 is used through an optical parametric oscillator to generate 
a squeezed vacuum state \Osq}, which is sent through BS-2 whose 
output is given by Eq. Jsl. The ^/-polarized and //-polarized output 
modes of BS-1 and BS-2, respectively, are sent through PBS-2 and 
PBS-3 to homodyne detection processes D-1 and D-2. The active 
retarders (or modulators) AR allow that only one of the V oi H po- 
larizations to proceed to D-1 and D-2 conditioned to photon-addition 
process in PA and on external control. The resulting state coming out 
from reconstructions in D-1 and D-2 is given by Eq. iHl. 



where 

|V„,o)=C/(V4)|n,0) = 



2"/2 ^ 
fe=0 



1/2 



\k,n-k), (3) 



are the eigenvectors of the density matrix p. {|V'n,o)} - with 
ri = 0,l,...,oo- form an orthonormal set, since the unitary 
U preserves inner product, and are also permutationally in- 
variant, i.e., given H = J2r,j \hj){j,^l nlV-n.o) = l^n,o>- 
This implies that Up = pH = p. The state (El alone cannot 
generate PPT-bound entangled states, as shown in reference 
II22I . We thus consider a more general class of states, given 
by 

p' = xp + {i-x)\n){n\, (4) 

corresponding to the mixture of state p in (J2]i with the state 

\n) = {i-\.f)y^±. h + 'l^^ .''\^,,,), (5) 

k=Q V '"^ 



whose preparation is depicted in Fig. 1 as we now explain. 
State (J5]l is generated by passing a //-polarized one-mode 
squeezed vacuum state. 



10. 




(6) 



generated in an optical parametric oscillator, through the 
beam-splitter BS-2. Here uj = (,^/|^|) tanh(|^|/2) is the 
squeezing parameter and thus < |a;| < 1. We impose for 
BS-2 that the parameter in ([T]| satisfy ^ 7r/4, i.e, not a 50 : 
50 beam-splitter. For simplicity, we take < 6 < tt/A. Thus 
we have \n) = U{0)\Osq), and writing |(/)„^o> = U{e)\n,0), 
we arrive at expression (Bl, with 
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Now to prepare the mixture in (HI, the two independent 
preparations are recombined in the polarizing beam-splitters 
PBS-2 and PBS-3 and proceed for homodyne detection on 
D-1 and D-2. The active retarders (AR) are externally con- 
trolled by the detection of a photon in the photo-addition pro- 
cess plus a external control to generate the full range of A in 
(pi. Whenever this photon is detected in the idler mode gen- 
erated in the parametric dow -conversion in BBO-1, indicating 
that a photon has been added to the signal mode, the AR al- 
low only the V^-polarized photon-added thermal component to 
leave to the detectors. When no photon is detected only the H- 
polarized two-mode squeezed vacuum is allowed to proceed 
to the detectors. By neglecting some of those detections, the 
ARs control the fraction A of polarization of the field incident 
at the photodetectors for repeated experiments. The idea be- 
hind this preparation is to use the fact that the beam-splitter is 
a classicality-preserving device and converts a (non-)classical 
state into a (in)separable one. The mixing in (J4l is thus a mix- 
ture of two nonclassical states and since the mixing parameter 
A is controllable by the experimentalist through the ARs, it is 
possible to prepare an entangled state whose partial transposi- 
tion is still positive. As we will see a crucial element here was 
the elimination of the vacuum component from the thermal 
state by the photon-addition process. 

Photon added and shifted thermal states. For the mixture 
present in Eq. (HI to allow a genuinely bound entangled state 
it is necessary to first change the Gaussian character of the V- 
polarized thermal state at mode A. A simple procedure to de- 
Gaussify the thermal state is the photon-addition as described 
in 1 29 1, which assumes that when the thermal state is fed as 
a signal into a parametric amplifier, the output signal state is 
conditionally prepared every time that a single photon is de- 
tected in the coiTelated idler mode. The simple assumption 
here is that the action of the conditioned parametric amplifi- 
cation is given up to first order in the coupling g between idler 
and signal by [1 + {ga\a\ — g*asai)\ where aV^ are bosonic 
creation operators acting on the signal (idler) modes. This 
results in a photon-added thermal state, which although pos- 
sessing the character needed (absence of the vacuum state). 



has failed to produce a PPT state. However if instead one 
is able to implement a saturated photon-addition |35| in the 
sense that the action of the creation and annihilation opera- 
tors in the signal is replaced by E^ = al{alas + 1)~^/^ and 
E~ = (ajfls + l)^^/^as respectively, the resulting state con- 
ditioned to the detection of one photon in the idler is given by 
the shifted-thermal state 
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where n is the mean thermal photon number These states 
were first considered by Lee in reference ll36l . in an analy- 
sis of the scheme depicted in |34|, where laser cooling sig- 
nificantly changes the emission of radiation of a micromaser 
The distribution dSll arises in 1 361 when the parameters of the 
micromaser cavity fulfill an ideal requirement. The shifting 
operation is also present in the related proposal of [37 1. Re- 
markably, a scheme to prepare ([8]l deterministically has been 
recently proposed fJSl, making the unconditional preparation 
of Q an achievable goal. The shifted thermal state allows the 
state (HI) to have genuine bound entanglement as we discuss in 
what follows. 

Entanglement detection. We now give a set of general rules 
for determining an arbitrary PPT bound entangled state. 



Observation 1 State Q is entangled for any < |a;| < 1. 



Proof: We use the following theorem, known as the Range 
Criterion 11171 : 



Theorem 1 For every separable state, there exists a set of 
product vectors {|a,6)} which spans the range of its density 
matrix, such that {|a, 6*)} spans its partial transposition. 



A violation of at least one of the conditions of this criterion 
implies entanglement, thus we prove that the range of p' in 
Q does not contain a single product vector. The range of p' is 
spanned by {|'0n,o)} - with n = 0, 1, . . . , oo - and \Q) as given 
in Eqs. (J3]l and (J5]l, respectively. We show now that assuming 
a product vector in the range of p' leads to a contradiction. 

Thus, let us assume that there exist complex numbers ruk, 
/c = 0,1,2,... such that 
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|f^)+^mfc|V'fc,o) = \vi)®\v2) 
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k=l 



write explicitly some important terms in left-hand side of ([9]l: 

mo|0,0)+mi(|0,l) + |l,0)) 
+ ((w/V2)mo cos^ e + ma) |0, 2) 

+ ((a;/\/2)mosin6'cos6' + \/2m2) |1,1) 

+ ((a;/V2)mo sin^ 9 + ms) |2, 0) 

+ ma (|0, 3) + Tail, 2) + ^312, 1) + |3, 0)^ 

+ ('moCj2(3/8)cos6lsin'^6l + m4V4) |1,3) 

+ fmow^ (3/8) 008^61 sin 61 + m4V4) |3, 1) 



We disregard normalization factors, since one can always in- 
corporate these factors in the values mi. The first terms in the 
right-hand side of (|9]l are 

ao/3o|0,0)+ao/3i|0,l)+ai/3o|l,0) 
+ ao/?2|0,2)+ai/3i|l,l>+a2/3o|2,0) 
+ ao/33|0,3) + ai/32|l, 2) + a2/3i|2, 1) + a3/?o|3,0) 
+ ai/33|l,3)+a3/3i|3,l) 

+ ... 

Let us first assume m^ — 0, this implies ao/?o — 0, by the lin- 
ear independence of the set {|j, j)} - with i, j = 0, 1, . . . , oo. 
Hence either ao = or /?o = 0; it is straightforward that 
either case would imply rrii — for all i, clearly being a 
contradiction. Let us then consider the case mo 7^ 0. With- 
out loss of generality, we assume ao = /^o = 1; from ex- 
pression (J9]), one gets mo — 1 and since (li/'n.o)} is per- 
mutationally invariant, for all n = 0, 1, . . . , cx), we have for 
the first odd terms mi = ai = /3i, and 7713 = a^ = f3^. 
So ai/33 = a3/3i. However, we have for the correspond- 
ing terms ai/Ss = mou;^(3/8) cos^sin'^ + m4-\/4, and 
Qg/?! = mow^(3/8)cos^6'sin0 + m4V4. Since 9 ^ 7r/4, 



we conclude that a^iPi 7^ ai/Ss, a contradiction [261. Thus, 
there is no product vector in the range of p' , implying Obser- 
vation 1 is true, QED. 

As pointed out in |19|, it is possible to construct entangle- 
ment witnesses for these states by applying the optimization 
methods of references [23 1 to projectors over the range of p' . 
Let us now check the conditions to obtain a PPT state. The 
following proposition relates the photocounting distribution 
of the input state pi to the PPT property of the output state 

Proposition 1 State pi is PPT if the following infinite- 
dimensional Hankel matrix l\24]l . 
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where |i;i) = X)fe=o'^fel^)' 1^2) = '^k'=oh'W)- Let us is positive semidefinite. 



Details of the proof are left for the Appendices A and B. A 
crucial point is that under a special ordering of the basis of the 
total Hilbert space Ha 'S) T-Lb, the partially transposed matrix 
of Q assumes a block-diagonal structure. 
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where each block can be decomposed as a Hadamard product 
of two matrices Mi — Aio B^, being both 

/ P\t\ P\i\ + 1 P\i\+2 ••• \ 
P\i\+1 P\i\+2 P\i\+3 ■■■ 
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and Bi positive matrices described in Appendix B. Since the 
Hadamard product preserves positivity and Bi is positive, p^^ 
will be positive semidefinite if Ai is positive semidefinite. But 
through Sylvester's criterion [301 one sees that this condition 
can be reduced to positive semidefiniteness of Aq. Thus the 
positivity under partial transposition can be checked by a hi- 
erarchical sequence of photocounting probabilities matrices. 
Indeed, a similar connection is already present in reference 
II22I with regard to the Stieltjes moment problem, in terms 
of nonclassicality exhibited by states Negative under Partial 
Transposition (NPT). An advantage of our approach is the 
direct-sum structure in Eq. (24i, which allows us to handle 



the construction of PPT states in a simple way. Similar direct- 
sum decompositions of the partially transposed density matrix 
can be found in ||20l I2TI and bring a great deal of simplifica- 
tion when dealing with PPT bound entanglement. 

The idea to obtain a PPT (HI is that the only block (111 for 
the shifted-thermal state that is negative under partial transpo- 
sition is the block AIq, due to po — 0. Hence, with a suitable 
small squeezing parameter w, the vacuum amplitude of the 
squeezed state, i.e |(00|risq)p = ^1 — |ajp will be suitably 
big in such a way to replace the null vacuum amplitude of 
p, making the mixture p' positive under partial transposition. 
The other terms are 0(|a;|), being as small as one desires, rep- 
resenting a small perturbation under control. We obtain thus: 



Although in practice, for continuous variable's regime it is an- 
alytically and numerically hard to determine exactly |27| the 
spectrum of p', it is possible to obtain upper bounds for the 
values of |w|, which guarantee a positive partial transposition. 
For the example considered, a conservative upper bound of 
\uj\ < lO^'^ has proved to be sufficient, allowing the gener- 
ation of a bound entangled state in continuous variables. It 
will remain an open question whether we have found or not 
a generic continuous variable bipartite bound entangled state, 
i.e., a bound entangled state with infinite Schmidt rank ||T9l . 

We have presented a procedure to unconditionally pre- 
pare bipartite bound entangled states in continuous variable's 
regime. The approach assumed was to deGaussify a ther- 
mal field by a photon addition process and to mix it with a 
squeezed vacuum state. Several minor results were developed 
in order to achieve this goal. Particularly the links with Han- 
kel operator theory and the direct-sum structure of the par- 
tially transposed state allowed us to give bounds on the pa- 
rameters that enable PPT bound entangled states to be pro- 
duced. There are few examples of such states in continuous 
variables and thus the novel class (Hli is interesting in its own, 
even if it were impossible to envisage a scheme to generate 
it experimentally. Our work shows that this is not the case 
and that the preparation in practice is feasible, opening new 
possibilities in quantum information processing protocols, as 
well as in the theory of quantum entanglement. We focused 
on thermal fields, due to thek implementation simplicity in 
the laboratory and also due to their use in fundamental ex- 
periments as |29|, but we stress that similar results could in 
principle be achieved with different classical photocounting 
distributions. Our approach was to "kill" the vaccum contri- 
bution of the thermal state by performing a photo-addition and 
then replacing it with the vaccum term of the squeezed state, 
which is superposed with other terms. If one is able to pro- 
duce a one-mode state satisfying Proposition 1 and then could 
entangle its vaccum term with other convenient state, we be- 
lieve it is possible to construct similar states to (Hh; this point 
is currently being investigated. 



Observation 2 There exists states Q which are PPT. By Ob- 
servation 1 these states are bound entangled. 

It is intuitive that Observation 2 is true given that the eigen- 
values of a matrix are continuous functions of its entries and 
thus a small perturbation of these values will not change the 
eigenvalues significantly. A constructive proof of Observation 
2 is left for the Appendix C, where we consider the example 
of such states with a balanced mixture A — 1/2 and a shifted- 
thermal state n = 1, that is 



p'^\{p + \n){n\), 



with 
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Appendix A: Block structure of output 

In terms of Fock basis, the state (|2| reads 

oo n 

P^T.P- E c^^,l\k,n-k){k',n-k'\, (15) 

n=0 fe,fe'=0 



where we introduced the symbols Cj,"^, = 2^v(fc)(fc')- ^® 
proceed now in order to find the structure of the partial trans- 
position of the density matrix p, which is obtained performing 
the operation of transposition in only one of the subsystems. 
From expression (J2]) and using permutation invariance. 



oo n 



P = T.P-^ E ci:'l\k,n~k){n~k',k'\, (16) 

n=0 fc,fc'=0 

hence the partial transposed matrix of the second mode in 
Fock basis is given by 

oo n 

P^^=Y.P^ E c'^kl\k,k'){n-k',n-k\. (17) 

n=0 k,k'=0 



An arbitrary element of this matrix is 

n 
(a, 6|/9^^|C, d) = 2jP« z2 '^k,k'^a.,kSb,k'5c,n-k'Sd,n-k, 



and we have then (a, 6|p-^^ \c,d) — unless a+d = c+b = n, 
or, equivalently, a — b = c — d. Taking this rule into account, 
we choose a special ordering of the basis of the Hilbert space 
H = Ha <8) Hb so that the matrix p^^ has a special block 
structure in this ordering. Let us first define the following 
sets: 

Bo = {\jj)}T=o, (18) 

B+,^{\i + k,k)}ZLi, (19) 

B^, = {\k,i + k)}^^,. (20) 

The notation should be clear: the elements in the set Bq are 
vectors \ab) which fulfill a — b — 0, while elements in set 
B±i are those which respect a — b = ±i. The union of these 
sets is precisely the basis of the total Hilbert space H, but 
now in a different ordering, i.e., we are ordering vectors \ab) 
according to their difference a — &. If we take as our ordered 
basis B = [J- Bi, the partially transposed matrix will show the 
block structure 



n^B 



/.. 








\ 




M_i 












Mo 












M+i 




V 








■ ■/ 



(21) 



and we write 



n^B 



= ^- 



(22) 



with each block Mi having the matricial representation 






/ PiCifi Pi+lCi^o Pi+2Ci,0 

(i+l) (i+2) (i+3) 

ft+lC-+l,l Pi+2Cl+i^{ Pi+3cl+ij 

Pi+2Ci^2,2 Pi+3'^1+2,2 Pi+i'^i+2,2 



\ 



\ 



■■/ 



where the notation X = means the matricial representation 
of X in the basis C. Thus, the state p will be PPT iff all the 
blocks Mi are positive semidefinite for all values i <E Z. 

To show the direct-sum decomposition in another way, we 
define T-Li = Span{Bi} {Span{.} amounts to the linear span 
of {.}) and we have the direct-sum decomposition of the total 
Hilbert space as 



H - H. 



(23) 



Appendix B: Proof of Proposition 1 

As shown in Appendix A, the partially transposed matrix of 
Q assumes a block-diagonal structure. 



0^^^ 



(24) 



Each block can be decomposed as a Hadamard product of two 
matrices M, = A, o Bi, with 



A, 



B, 



p\i\ Pm+i p\i\+2 

P\i\ + 1 P\i\+2 P\i\+3 

P\i\+2 P\i\+3 P\i\+4: 

^,0 ''4,0 ''i.O 

(i+l) (i+2) (i+3) 

^-1+1,1 '-i+l.l ^i+1.1 

{i+2) (i+3) (j+4) 

'-i+2,2 ''i+2,2 ^+2,2 



\ 



(25) 



V 



.•\ 



■■■/ 



(26) 



We will need the following theorem, known as the Schur 
Product Theorem |30|: 

Theorem 2 The Hadamard product of two positive semidefi- 
nite matrices is a positive semidefinite matrix. 

To prove Proposition 1, we prove first a lemma: 

Lemma 1 The matrices Bj are positive definite, for all j £ Z. 

(fe) _ 



Proof. First, starting from the definition c! 



2fc\/iJij7 ~ 2*^ Y 



^J 



i\j\(k-iy.{k-jy. 



, we express an arbitrary 



Bj as a Hadamard product of six matrices: 



B, 



C-i o Dj o Ej o e] o Fj o F^ , 



(27) 



with 



c,= 



D,= 



E,= 



( j! (i + l)! (j+2)! 
(j + 1)! (j + 2)! (j+3)! 
(j+2)! (j+3)! (j+4)! 

V ; ; ; 

/ 1/2^' 1/2^+1 1/2J+2 .. 
1/2^+1 1/2^+2 l/2J'+3 .. 
1/2J+2 1/2^+3 1/2J+4 .. 

V ; ; ; ■ 

/l/VO! 1/^ 1/VO! ...~ 
l/\/lI l/\/lI l/\/II ... 
1/V2! l/%/2l 1/V2! ... 



The subspaces Hi are invariant under the action of the oper- 
ator p^^ , i.e., this operator does not send vectors from Hi to K 
a different Hj. So, the operator p ^ should decompose as a 
direct-sum l.33i . 



/ 1/^ 



■■/ 



1/V7 



1//J 



l/x /(j+^ l/x /(J+^ l/v /(J+I)! 
l/v/(jT2)! 1/7(7+2)! 1/7(7+2)! 



V 



■•/ 



The matrices Dj, Ej and Fj are rank-1 matrices, so are ob- 
viously positive. The Hankel matrix Cq is positive definite, 
since its leading principal minors of order k have determinant 
n*^^Q(i!)^. A similar argument holds for an arbitrary Cj and 
thus they are all positive definite. Another way to prove this 
is observing that the following sequence satisfies the Stieltjes 
moment problem 11311 : 



/„ = n!=/ a;"e~=" dec, n = 0,1,2,... (28) 
Jo 

So, Co is positive definite. Now, we see that 



(29) 



Jo 

/"OO 

gi = {j + iy.= x^{x^e-'')dx, (30) 

Jo 

/•OO 

92 = {j + 2)\= / x^ix^e-^dx, (31) 

Jo 

: (32) 

defines a sequence that also satisfies the Stieltjes moment 
problem. From Theorem 1, the Hadamard product of positive- 
semidefinite matrices is a positive-semidefinite matrix. Also, 
the Hadamard product of a positive-definite matrix - Cj in 
our case - with rank-1 matrices - Dj, Ej and Fj - is positive- 
definite ESJI . So Bj is a positive-definite matrix, QED. 

The proof of Proposition 1 is now straightforward. Since all 
Bi are positive definite, to have all Mi positive semidefinite 
we must have all Ai positive semidefinite. But Aq positive 
semidefinite implies, by Sylvester's Criterion, that all other Ai 
are positive semidefinite, since they are principal submatrices 
of Aq. So, all Mi = AiO Bi are positive semidefinite if Aq is 
positive semidefinite and by the block structure of p the state 
is PPT if Aq is positive semidefinite, QED. 



Appendix C: Proof of Observation 2 

For simplicity, we will construct a example of a PPT (|4| 
with A — 1/2 and p being the output of (l8| with n = 1. Thus, 
the state we are considering is 



with 



p' = -ip+\n){n\) 



OO / 1 \ ^ 

XI ( 9 ) IV'"o)(V'«o| 



(33) 



(34) 



We first observe that any shifted thermal state p generated by 
(|8]l is locally equivalent to ( [34| above. Define the following 
invertible operation; 

T\n) = in + iy/'f^Y'\n) (35) 



Then we have that T (g) TpT^ (g) T^ is equal to the matrix 
( 34 1. Local invertible operations do not affect positivity of the 
partial transposed matrix, thus this special case is broad in this 
sense (for the case of A = 1/2). 

We know that ( 34 1 is NPT, since po — 0, implying that pi 



is nonclassical and thus NPT, by the criterion of 1221 . Also, 
from Appendix B, we know that 



-e 



M,, 



(36) 



where Mi — AioBi. For i 7^ 0, the blocks. Mi are all positive 
definite, since the corresponding matrices Ai are positive def- 
inite. When we consider the new partially transposed matrix 
for ( |33| ), we have 

p'^-=(l/2)(p^« + |0)(r!r-). (37) 

We will now consider that |a;| is sufficiently small that we can 
neglect terms 0(|aj^ |); we will discuss more on this point fur- 
ther in this Appendix. Thus, in this approximation we can say 

that effectively we have I r2)(f] I "^^ « (^1 - |a;|2)(|G0)(00| + 
'^l02o)(OO| +w|00)((/)2ol)- We can rewrite Eq. ([37]i as 



„'T„ 



(1/2) 



e^4 



p 



(38) 



and A/^ = Mq 



where Af/ = M„ for i ^ 
v/l-|^P|00)(00|, while P = c^(|02o)(OO| + |OO)(02ol)^" 
represents a perturbation with magnitude totally dependent on 
the value of |w|. 

It is straightforward that for values -^/l — |wp > 1/2, 
which means |a;| < ■y/3/4, the block A/q becomes positive- 
definite. We must know how small |a;| must be so that p'^'^ 
remains positive. 

Remark: Since the eigenvalues of a matrix are continuous 
functions of its elements, we could already stop the demon- 
stration at this point: slight variations of uj would not affect 
the eigenvalues of the positive blocks M- significantly and 
consequently its positivity. However, the constructive demon- 
stration given here has the advantage of giving an estimate 
on the order of magnitude of the value cu, which is relevant 
experimentally. 

We need the following theorem (Theorem 6.1.1 from ll30l ). 
known as the Gerschgorin Disc Theorem: 

Theorem 3 Given a n x n square matrix M = [rriij], let 

n-l 

R^iM)= Y. I"^*jl' * = 0,1,..., 77-1, (39) 

3=0, j^i 

denote the deleted absolute row sums of M. Then all the 
eigenvalues of M are located in the union ofn discs 

n-l 

\J{zeC:\z-a,,\< R,(M)} = G{M). (40) 

i=0 

Furthermore, if a union ofk of these n discs form a connected 
region that is disjoint from all the remaining n— k discs, then 
there are precisely k eigenvalues of M in this region. 
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The region G{M) is called the Gerschgorin region of M and 
the individual discs in G{M) are called Gerschgorin discs; the 
boundaries of these discs are called the Gerschgorin circles. A 
similar result holds for the collum-sums (Corollary 6.1.3 from 
|30|), but since we are dealing with Hermitean matrices it will 
not affect the results. We need also the following refinement 
(Corollary 6.1.6 from |30|): 

Corollary 1 Let D = diag{d^, di, . . .}, with di positive real 
numbers for alii — 0,1, .. .. Then all eigenvalues of M lie in 
the region 



[J {z G C : |z - a^i\ < 



n-l 



i=0 



1 

di ^-^ 



dj\mij\} 



G{D-^MD). 



Moreover, the spectrum of M is precisely the set 
Hz, G{D-^MD). 

Since the blocks M[ are all positive-definite, there exists a 
set of diS above which will bring all Gerschgorin disks to the 
positive segment of the real line. Also, there is a continuous of 
such di's and we conclude then that a slight change in the row 
sums Ri - which are constituted by the off-diagonal elements 
of the matrix - will not change the eigenvalues of a matrix. 



since this would correspond to a negligible deformation of the 
corresponding Gerschgorin region. Let us see how this applies 
to example ( [33] l. 

The perturbation matrix P affects only the first row and 



collum of each A/q, 


M^ti ■' 


ind A/±2- For the first fines of these 


blocks we have 


ko 


-Vi 








-a;2| < 2a; + iT!o 


(41) 




ki 


Pll 

2 ' 


<uj + Ri 


(42) 




\Z2 


4 ' 


<UJ + R2 


(43) 



By putting the sole value uj, instead of the actual terms 
ujcos^{9), ujsm^{9) and a;sin(6') cos^(6'), we are doing an 
overstimation of the perturbation. By the equations above, if 
1^1 << P2/4 < pi/2 < Y^r^H^' we will have that p''^^ 
will remain positive, since the associated Gerschgorin regions 
will be effectivelly unaffected. The value P2/4 = 1/16 has 
a order of magnitude of 10^^; thus we will impose for |a;| a 
conservative upeer bound of 10^"^. We can now justify the ne- 
glecting of terms O(w^): their rate of decrease is much faster 
than the rate of decrease of diagonal elements; also, we have 
not considered terms sin'''(0) cos^^*^(0), which would make 
this rate even faster. 



